Spiral order induced by distortion in a frustrated square-lattice antiferromagnet. 
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In a strongly frustrated square-lattice antiferromagnet with diagonal coupling J', for a — 
J/(2J') < 1, an incommensurate spiral state with propagation vector Q = (tt ± 5, tt ± 5) near 
(7r, tt) competes closely with the Neel collinear antiferromagnetic ground state. For classical Heisen- 
berg spins the energy of the spiral state can be lowered as it adapts to a distortion of the crystal 
lattice. As a result, a weak superstructural modulation such as exists in doped cuprates might 
stabilize an incommensurate spiral phase for some range of the parameter a close to 1. 

PACS numbers: 75.10.-b 75.25.-fz 75.50.-y 75.90.+W 



An interplay between small distortion of the crystal 
lattice and the magnetic properties of the material is 
currently a subject of intense research. One problem 
which supplies strong motivation for such studies is that 
of stripe order in the lightly doped high-T c cuprates 
La2_ 2; Sr a ;Cu04+j ) (LSCO) and in related nickelates 
These phases are always associated with a weak super- 
structural distortion of the original "stacked square lat- 
tice" structure of the un-doped parent material. Incom- 
mensurate magnetism in these compounds is usually in- 
terpreted in terms of a segregation of the doped charges 
into lines which separate the antiferromagnetic domains 
( "stripes" ) characteristic of the un-doped material. Al- 
though modulation of the crystal structure which is in- 
duced by charge-stripe segregation is often too small to 
be observed in experiment 1], it is clear that essential 
result of the stripe order for the spin system of cuprates 
is a periodic modulation of the exchange coupling in the 
Heisenberg spin Hamiltonian which describes their mag- 
netic properties [3j. So far, though, only the simplest "av- 
erage" consequence of the stripe superstructure, in the 
form of the effective weakening of exchange coupling in 
the direction perpendicular to the stripes, has been con- 
sidered g|. A similar problem, of an interplay between 
the spin order and the cooperative Jahn- Teller distortion 
accompanying the charge order, arises in the context of 
the charge-ordered phases in doped manganites || . 

Because the low-energy magnetic properties of layered 
LSCO cuprates are believed to be adequately described 
by the two-dimensional (2D) Heisenberg spin Hamilto- 
nian, this model has recently become a focus of intense 
research. Special attention was devoted to the frustrated 
square lattice, where in addition to the nearest-neighbor 
exchange interaction, J > 0, there is a diagonal coupling, 
J' > 0, such that a = 577 is close to 1. It was originally 
motivated by the predictions that non-Neel resonating 
valence bond states 0, and quantum-critical behavior 
H associated with the T — order-disorder phase tran- 
sitions which may occur in this case might be important 
for the physics of the superconductivity in cuprates. 

Despite RVB spin-liquid state and quantum critical- 
ity are strongly predicated upon the quantum nature of 



the spins (S=l/2 in cuprates), a semiclassical spin- wave 
theory appears to provide a surprisingly good guidance 
to the behavior of the frustrated s qua re-lattice antiferro- 
magnet (FSLA) H HE1 HHlll Hi. Perha P s < this 
is because the phenomenon of frustration mainly rests 
on the ground state degeneracy which exists for classi- 
cal, as well as for quantum spins. In fact, existence of 
the spin-liquid phase possibly related to the RVB state 
in the FSLA for the range of the parameter a around 
a = 1 was first conjectured in Ref. |3 on the basis of the 
conventional spin-wave calculation to the order 1/S. This 
suggestion was then supported by the field-theory meth- 
ods (l6 j. n umerical calculations [ToL fill ll2T ] and other 
studies |l3l [lj, |l5j . It was established that a disordered 
phase, whose nature is still controversial, is realized for 
0-8 ;$ ct < 1.1. Although these studies were essentially 
aimed at understanding the physics of doped LSCO and 
related materials, the lattice modulation was generally 
ignored. One reason for this is that, traditionally, the lat- 
tice distortion in a spin system is treated by switching to 
a larger unit cell, with multiple different spin sites. This 
approach is not viable for the long-periodic superstruc- 
tures, and is not possible for the charge-ordered states 
with incommensurate modulation. 

In contrast with the previous studies, present paper 
addresses the consequences of the superstructural lattice 
distortion for the ground state of the 2D Heisenberg spin 
Hamiltonian with classical spins, ie essentially presents 
an "unrealistic" mean-field (MF) treatment of the realis- 
tic spin model. Although MF results are subject to sig- 
nificant quantum corrections, especially for small spins 
S< 1, they nevertheless provide useful guidance about 
the hierarchy of the competing ground states (GS) in the 
system. In fact, the MF ground state very often survives 
account for quantum and thermal fluctuations, as it does 
for the un- frustrated 2D antiferromagnet. 

Main finding of this paper is that a weak superstruc- 
tural modulation of the crystal lattice in the FSLA 
may stabilize an incommensurate spin-spiral ground state 
with the propagation vector Q — (tt ± 5, tt ± S) close to 
(tt, tt) for a < 1. Although in the absence of a struc- 
tural modulation the energy of the spiral states is higher 
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than that of the collinear Neel states illustrated in Fig. 
1(a) (except for a = 1), they are in close competition 
for a near 1. While spin spiral is usually ignored in the 
analysis of the possible phases in FSLA, in presence of a 
superlattice modulation it might actually win the com- 
petition for some range of a < 1. Here this is shown 
explicitly on the mean field level, by treating the effect 
of a small but otherwise quite arbitrary lattice distor- 
tion, as a perturbation in the microscopic classical-spin 
Heisenberg Hamiltonian. 

Consider a system of N equivalent spins on a square 
lattice, Fig. 1 (a), coupled by Heisenberg exchange in- 
teraction, H = J2i j Jij (SiSj)- While only coupling be- 
tween the nearest neighbors along the side (J) and along 
diagonal (J') will be of interest in this paper, here = 
Jji parameterize a general exchange coupling between 
the spins at arbitrary lattice sites i and j. In the absence 
of a distortion the MF classical ground state is a planar 
transverse spin spiral, Sj — (S cos(Qrj), Ssin(Qrj), 0), 
fl?l Hif . The ordering wave vector Q corresponds to the 
minimum of the lattice Fourier transform of the exchange 
interaction, Jq = ^2 r . . Jij ; exp(— igr^), = Tj — r». 
This GS is obtained by finding the minimum-energy con- 
figuration for the Heisenberg Hamiltonian with classical 
spins under the constraint that S 2 — S 2 for all sites j. In 
general case, spontaneous symmetry breaking is defined 
by the two mutually perpendicular spin vectors which de- 
termine the polarization of the spiral, ie by the Fourrier 
transform of the lattice spin distribution, Sq = S' +iS" . 
For collinear situations, such as ferro- or antiferromagnet, 
corresponding to Q — and, eg, Q = (w, n), respectively, 
only a single vector is needed for the order parameter. 

A slight distortion of the crystal structure which is 
characterized by the appearance of the additional, weak 
supperlattice Bragg reflections at wavevectors ±Q C , cor- 
responds to a small harmonic modulation of the ionic 
positions, (tv,)' = rj + eiCos(Q c r 3 ) + e 2 sm(Q c r., ). In 
most general case, this results in a harmonic modula- 
tion of the exchange coupling. It has either the same 
wavevector Q c , if it appears as a first-order correction to 

Jij in small parameter e ~ f^r J < 1, or the wavevec- 
tor 2Q C , if it appears only in the second order, ~ e 2 , [l9| . 
There is also a second-order correction to the bond en- 
ergy, — + 5 Jij. The spin Hamiltonian becomes, 



a) 



H = 



i,3 



SiSj. (1) 



where the tildes were omitted, and the complex j {j — 

Jij + x 3ij was introduced. While without distortion Jy 
would satisfy all symmetries of the lattice, exchange con- 
stants in Eq. 0) possess only those symmetries of the 
un-distorted lattice which preserve Q c and the distortion 
polarizations 61,62 (this includes all translations). 

The modulated-exchange terms allow um-klapp pro- 
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FIG. 1: (a) Frustrated square lattice with diagonal coupling 
J/(2J') < 1; continuous degeneracy of the MF ground state 
corresponds to an arbitrary angle between the two antiferro- 
magnetic sublattices. (b) "Stripes" on the square lattice with 
diagonal modulation, Q c — ( „ j ^r)j n = 4 case is shown. 

cesses which couple Sq and Sg±Q in the spin Hamilto- 
nian, and couple these Fourrier-components in the equa- 
tions expressing the conditional minimum of the classical 
exchange energy. As a result, additional Fourrier har- 
monics, at wavevectors Q + nQ c , n = ±1, ±2, appear 
in the GS spin structure. It has the form of expansion, 



Sq = J2S 



Q+nQ c °q,Q+nQ c + b Q+nQ° q,Q+nQ c 



where Sq ~ 0(e'™' ). This corresponds to a bunched 
spiral 0, llSlIT^ . Based on very general exchange sym- 
metry arguments 0, 0, H3] , in the absence of any ad- 
ditional symmetry breaking, the perturbing terms have 
to be proportional to the non-perturbed order parame- 
ter. As a result, the leading new Fourrier-components, 



'Q+Qc 



y Q-Q c 



X±(Q c )u 2 Qc S-z 
X±(Q c )"qS- 2 



S Q +0(e 3 ), (2) 



S Q +0(e 3 ). (3) 



Here u q = S^2(J q - J Q )(J q+ Q + J q q - 2J Q ) is 
the spin-wave spectrum in the initial, non-distorted, 
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single-Q exchange spiral, X±(q) = %(Jq — Jq) 
is its transverse (perpendicular to the spin plane), 
g-dependent staggered static spin susceptibility, and 
jq = J2vi 3ij ex P(~ l 1 r ij) = j-q 1S a lattice Fourrier- 
transform of the modulated exchange term. Neglecting 
the 0(e 4 ) terms, the corrected ground state energy is ob- 
tained in the second order of the perturbation theory, 



cos2Q'). Upon account for distortion the GS energy is, 



K 



GS 



N 



= JqS 2 - 



3Q-lQ«-3QHQc 



S 4 



X±(Q>^ c 



(4) 



and, correspondingly, is, in general, lowered by the ex- 
change modulation. This occurs as a result of the appro- 
priate adjustment (bunching) of the initial single-Q spi- 
ral spin structure, through appearance of the additional 
Fourrier-harmonics, Sq ± q . In addition, the pitch of 
the primary spiral component, Sq, may also change, 
Q — > Q, because the spiral propagation vector, Q, is now 
defined by the minimum of the corrected energy, Eq. 

A singular situation occurs when the lattice modula- 
tion has the wavevector Q c which is near the dispersion 
soft spot of the initial spiral, eg close to its Goldstone 
mode. In this case corrections © - diverge, and 
the perturbation approach fails, highlighting the sensi- 
tivity of spin system to such modulations. In frustrated 
spin systems, entire soft regions, such as the lines of 
soft modes, often appear due to the accidental cancel- 
lation of the interactions. As a result, such systems must 
be extremely sensitive to structural distortions. On the 
other hand, in many important cases, such as a nearest- 
neighbor non- frustrated antiferromagnet, Q is a special 
symmetry point of Jq {]q ~ £ Jq)i an d the correc- 
tion term vanishes. Therefore, simple structures, such 
as collinear antiferromagnets, are, in general, not sensi- 
tive to small lattice modulations. In what follows, the 
singular cases will be excluded from the consideration. 

The results obtained above can now be applied to an- 
alyze the effect of lattice modulation in a square-lattice 
antiferromagnet, which may be of direct relevance for the 
charge-ordered phases in doped LSCO cuprates and re- 
lated perovskites. For definitiveness, consider the case 
of n-periodic diagonal modulation with Q c = 
where n = 2,3,4,..., illustrated in Fig. 1(b) (in LSCO 
the most stable superstructure occurs for n = 8, Q). 
Without frustration, J„ = eJq, and, upon switching to 

■, the problem is fac- 



_ (a 1 -a 2 



Q = q iOi+Oal and Qi _ , ( T 

torized and corrections are essentially the same as for ID 
chain. There is no change of the global minimum of clas- 
sical spin energy, so the nearest-neighbor antiferromag- 
netism is stable with respect to the bond modulation. 

In the frustrated case, Jq = 4 J cos Q cos Q' + 
2J'(cos2Q + cos2Q'), and, if both side and diagonal 
bonds are modulated, jq — 4 J cos Q cos Q'+2j (cos2Q+ 



E 



sin Q 



j cos Q' + 2j' cos 2 cos Q 



NS 2 Q J cos Q cos Q' + 2 J' cos 2 cos 2Q 



(5) 



In the absence of bond modulation, the ground state 
is determined by the hierarchy of the local minima 
of Jq, which only depends on a — ^jj. For weak 
frustration, a > 1, the global minimum is that with 
sinQ = sinQ' = 0. It corresponds to the conventional, 
collinear Neel antiferromagnetic order with a single prop- 
agation vector Q = (jr, tt), and the ground state energy 



-4J5 2 (1 — Although there are four 



equivalent Q-points in the Brillouin zone (BZ), (±tt, ±7r), 
(=p7r, ±7r), which restore the lattice C4 rotational symme- 
try, they are related through addition of the appropriate 
reciprocal lattice vectors t, so there is no true GS de- 
generacy in Q-space. The only degeneracy is the GS 
rotational symmetry in spin space, which is left of the 
0(3) symmetry of the Heisenberg spin Hamiltonian. 

For strong frustration, a < 1, there are two non- 
equivalent lowest-energy minima of Jq, they satisfy 
cosQ = cosQ' = and have the GS energy j?E( w ,o) = 
—4J'S 2 = —AJS 2 ^. They correspond to two pairs of 
equivalent Q-points in the BZ, (±7r, 0) and (0, ±7r), which 
represent the antiferromagnetic order propagating along 
the X and Y axis, respectively. This double degeneracy 
in Q-space can be used to construct a continuum of states 
which are the linear combinations of the above two. This 
continuous GS degeneracy is usually described in terms of 
two decoupled antiferromagnetic sublattices based on the 
diagonals of the original square lattice, which is transpar- 
ent for J' ;§> J. Each sublattice has an antiferromagnetic 
order, but there may be an arbitrary angle between the 
two, because the mean field from one sublattice cancels 
on the sites of the other, Fig. la. This continuous degen- 
eracy is lifted by zero-point or thermal spin fluctuations 
which prefer collinear arrangements of the two sublat- 
tices in the GS. This phenomenon is known as "order 
from disorder" jHHl. 

Although it is not the focus of this paper, an interesting 
situation occurs for a = 1, when, on the MF level, there 
is also a continuous GS degeneracy in the Q-space. The 
minimum condition for Jq becomes cosQ = cosQ', and 
is satisfied for any spiral with the propagation vector Q 
that belongs to the square with the vertices at (±7r, ±7r), 
(=F7r, ±7r). They all have the same energy, jjE a=1 = 
—2JS 2 = —AJ'S 2 . This continuous Q-space degeneracy 
is at the origin of the spin-liquid phase conjectured in 
FSLA for a close to a = 1 0, E El E El Q E 13 • 

What is important here, is that the spiral states with 
Q w (n, it) are in close competition with the collinear 
states for a < 1. In particular, the spiral with the prop- 
agation vector defined from Q' = Q, cosQ = —577, ie 
Q = (cos _1 (— yjt), cos~ 1 (— 577)), is a local minimum of 
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Jq along the diagonal, (q,q), direction, whose energy in 
the absence of modulation is jjEq = —2aJS 2 . Except 
for a = 1, though, the energy of this local extremum 
(and of all other spiral states) is higher than that for the 
decoupled antiferromagnetic sublattices, E(~m, and for 
this reason they are usually ignored. However, it is clear 
from the Eq. (J5J that, while the energy of the antifer- 
romagnetic states is insensitive to the bond modulation, 
the energy of the spiral state can be lowered as it adapts 
to the lattice distortion! Therefore, at least on the MF 
level, a spiral may become the lowest energy state (ie 
the ground state) for some range of the parameter a in 
the vicinity of 1 (whose width is ~ 0(e 2 )). For a long- 
periodic modulation, Q c <C 1, and for j = 0, it is easy to 
find that the spiral phase is stable for 1— \jj J| 2 < a < 1. 
The principal propagation vector of the spiral is obtained 
by minimizing Eq. JSJ). 

While it would be interesting to study the modulated- 
exchange Hamiltonian JQ) for quantum spins and for the 
arbitrary values of \jij/J%j\, this is a formidable task 
which is beyond the scope of this paper. Here a per- 
turbative scheme is used to find the mean field ground 
state. It is valid for classical spins, S > 1, and for 
small exchange modulation, Ij^/Jij] ~ e <C 1- Never- 
theless, it provides an important insight into behavior 
of the frustrated square-lattice antiferromagnet. A find- 
ing that (by selecting the spiral order) exchange mod- 
ulation effectively destabilizes collinear Neel states pre- 
ferred by the fluctuations clearly supports the instabil- 
ity of the frustrated square-lattice antiferromagnet with 
J/(2J') close to 1 with respect to the bond-modulated 

states, mm [13. 

The essential results of this paper are summarized by 
Eqs. J5J - Q- The main finding is that the energy of 
the equal-spin transverse spiral state can be lowered by 
the exchange modulation in the Heisenberg spin Hamil- 
tonian. This happens as spiral adapts to the modu- 
lation through appearance of the additional Fourrier- 
harmonics, Sq +k q , n = ±1,±2, ... (bunching). As a 
result, in frustrated square-lattice antiferromagnet with 
diagonal coupling J', such, that a — J/ (2 J') is close to 1, 
lattice modulation may open a region of stability of the 
incommensurate spiral phase. This "order by distortion" 
phenomenon competes with "order by disorder" , which 
prefers collinear arrangements of two antiferromagnetic 
sublattices. Incommensurate spiral phase with the prop- 
agation vector Q = (ir ± S, tt ± S) close to (ir, ir) wins for 
the range 0(e 2 ) of a around a = 1. 

The arguments presented here provide plausible ex- 
planation for the incommensurate spin-ordered phases, 
which are among the most interesting and puzzling fea- 
tures observed in the doped perovskites, and may also be 



of direct relevance for the doped LSCO materials. For the 
Heisenberg spin Hamiltonian on square lattice in the ab- 
sence of distortion, one needs at least a third-neighbor 
coupling in order to stabilize the MF spiral ground state. 
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